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0\ ■ Abstract 



The structure of the QFT expansion is studied in the framework of a new "Invariant analytic" version of 
the perturbative QCD. Here, an invariant (running) coupling a(Q 2 /A 2 ) — /3ia s (Q 2 )/47r , is transformed 
into a "Q 2 -analytized" invariant coupling a an (Q 2 / A 2 ) = A(x) , which, by constuction, is free of ghost 
singularities due to incorporating some nonperturbative structures. 

Meanwhile, the "analytized" perturbation expansion for an observable F, in contrast with the usual 
case, may contain specific functions A n {x) = [a n (x)] an , the "n-th power of a(x) analytized as a whole", 
instead of (A(x)) n . In other words, the pertubation series for F(x), due to analyticity imperative, may 
change its form turning into an asymptotic expansion a la Erdelyi over a nonpower set {A n {x)} . 

We analyse sets of functions {A n (x)} and discuss properties of non-power expansion arising with 
\^ ■ their relations to feeble loop and scheme dependence of observables. 

The issue of ambiguity of the invariant analytization procedure and of possible inconsistency of some 
, of its versions with the RG structure is also discussed. 

oo : 

1 Introduction 

^— > , 

In papers |JJ [| the procedure of Invariant Analytization of perturbative quantum chro- 
modynamics (pQCD) has been elaborated. It implements a combining of two ideas: the 
RG summation of leading UV logs with analyticity property in the Q 2 variable, imposed 
.> ; by spectral represent at ionP] of the Kallen[[|] - Lehmann[[|] type - see Eq.(|I|) below. This 

^ ■ combination was first proposed and devised || to get rid of the ghost pole in QED about 

forty years ago. 

The central notion of the RG improved QFT perturbation technique is an invariant 
charge or "invariant coupling" a(Q 2 /A 2 ) involved into the renorm-invariant expressions 
for observables. The QCD scale parameter A related to the ghost singularity position 
equals approximately to 300 MeV. For properly normalized QCD invariant (running) 
coupling a(Q 2 / A 2 ) = j3ia s (Q 2 )/An] (3i(n) = 11 — (2n/)/3 (defined in the space-like 
region), the analytization results in a specific transformation into a form a an (Q 2 / A 2 ) free 
of ghost singularities. Here, by construction, 

the analytic coupling is defined via the Kallen-Lehmann representation 



(a) da 

/ with the spectral density, pier) = 3 a(—a) , (1) 

7T J a + x 
o 



1 To appear in the April issue of Teor. Mat. Fizika v. 119 (1999) No.l 

2 This representation implements^] general properties of local QFT including the microscopic causality 
one. In this context we shal use the term "analyticity" as a synonym of "spectrality" and "causality" . 
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calculated on the basis of "initial" RG- summed perturbation expression a(x) for invariant 
coupling. 

Generally, a an (x) differs from its "original input" a(x) by nonperturbativeQ additive 
terms, which "subtract" unphysical singularities - see, e.g., Eqs. (fj) and (|6]). 

A detailed analysis revealed |2| that the analytic coupling a an (x) obeys several 
important properties. It turns out to be remarkably stable in the IR region, at Q < A, 
with respect to higher-loop contribution and to renormalization scheme dependence. Its 
IR limiting value a an (0) is universal in this sense. We review this subject quite shortly in 
Sections [H] — |2j. 

On the other hand, the Invariant Analytization of a physical amplitude F(Q, a) is 
not a strightforward procedure. A few different scenarios are possible. In papers [§. [jj, 
a particular version, the Analytic Perturbation Theory (APT), has been proposed and 
elaborated. Here, due to specific analytization ansatz, instead of the power perturbation 
series common for theoretical physics and QFT, an analytic amplitude ^(x) is presented 
in a form of an asymptotic expansion of a more general form, the expansion over an 
asymptotic set of functions A n (x) = [a n (x)] an , the "n-th power of a(x) analytized as a 
whole" . In the APT approach, the drastic reduction of loop and renormalization scheme 
sensitivity for several observables has been found - see Refs.p|-|10||. 

To understand the nature of the "APT's loop and scheme immunity" , in Section [3] we 
study properties of a nonpower asymptotic set {^4 n (x)} emerging from the APT recipe 
of analytization. 

In Section |], we analyse the structure of possible variants of analytization of expression 
for an observable and discuss the danger of inconsistency (more precisely — incompati- 
bility with the inner structure of RG) for some of them. 



2 RG solution and analytization 
2.1 One— Loop Analytization of a(x) 

At the one- loop level, the invariant coupling a(x) = 1/ lux suffers from a pole singularity 
at x = 1 incompatible with the spectral representation ([]]). Here, analytization consists 
of analytic continuation of the 1/ In a; expression into the negative x region and defining 
the spectral density via its imaginary part. The resulting spectral integral (|XJ) with p(a) = 

(in 2 cr + 7r 2 ) can be calculated explicitly 

1 11 

a i X ) = 1 a an(x) = A(x) = . (2) 

ma; in x x — 1 

The second term, precisely compensating the ghost pole, has a nonperturbative nature. 
At small a M = a(Q 2 = /i 2 ) it is not "visible" in Taylor series as it behaves like exp(— l/a M ) . 

3 On a deep connection between renormalization invariance plus causality and nonanalyticity in a, 
see our papers 0. 
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To see this clearly, one should return from A-parameterization to the one in terms of a M , 
the renormalized coupling constant, and /x 2 , the reference momentum squared. 

Note that a relation between A//i and a M in the course of analytization transformation 
changes. Instead of the usual expression A 2 = /z 2 exp(l/a M ), according to Eq.(fJ), we have 
the transcendental relation 

A 2 - ■ -- 1 1 1 



/(a)' In f(a) l-f(a) 

Here, at small a (as well as in (||]) at large x) , one can neglect the second, nonperturbative, 
term. Meanwhile, at a ~ 1 (x < 1) this term dominates, providing the IR fixed point at 
a = 1. 

The analytic coupling (|2]) is a monotonous function in the whole interval (0, oo) with 
the finite IR limit. The second term quickly diminishes as x — > oo: it contributes about 
5% at Q — 10A and only 1% at Q — 25A. The whole set of solutions @ with various A 
values, considered at the Q 2 scale, forms a bunch with a common limiting point a(0) = 1. 
This value, corresponding to a an (0) = 47r//?i Q, turns out to be universal. It does not 
change in the two- and three-loop approximation as well - see Refs. ffl, 0, [12[ . 



2.2 Two- and three-loop cases 

For the two-loop case, the invariant coupling has to be defined by the transcendental 
relation 

^- Mn ( 1 + 6^) =lnX; h= % ( = S atn/ = 3 ) (3) 

resulting from integration of the two-loop RG differential equation. 
The iterative procedure yields the explicit approximate solution 

used in our previous papers. 

The exact solution to Eq.fl3|) can also be expressed [JTTJ, [12 



a^\x) = — h - 1 + yv{x) ; W(x)=W_ 1 (z); z = -e-^ 1 (5) 

in terms of a special function W, the Lambert function (we use the notation of paper 
0) 

W{z)e w(z) =z, 

with an infinite number of branches W n (z). Most part of the physical region x > 1 
corresponds to the particular branch W-i(z) that is real and monotonous for z < — 1/e. 

4 The effective flavour number at residue of the pole, evidently is rif = 3. This gives a an (0) = 1.396. 



3 



The relation between z and x in Eq.([5|) is normalized to correspond with approximation 
(|) at a; > 1. The branch W(x) = W-\{z) is complex below the ghost singularity at 
x = 1 , z = — 1/e. 

The qualitative analysis of Eq.(|3|) or its exact solution shows that at x — 1 + \e\ the 
ghost singularity has a form of the square-root type branch point 



,( 2 ) 



1 



cr ' ( x 



2b{x - 1) 



— + 0{Vx~^T) 



3b 



that yields an unphysical cut between this point and the origin. 

To illustrate, proceeding from the two-loop solution (g), we define the analytic coupling 
by the spectral representation (|l|) with spectral density defined in terms of W, that is 
equivalent to subtraction of the "cut integral" related to the square-root singularity 

Note that iterative approximation (|]) has a slightly different structure of ghost singula- 
rities^]. 

Nethertheless, in Ref. [[Ll|] it has been demonstrated that the analytized iterative solu- 
tion is numerically very closef] to the analytized exact one, Eq.(||). As a practical result, 
this means that for the two-loop a an (x) one can use an expression in the form Eq.(Jj]) with 
spectral density 



R{L) = L + bln\l (l + lf) +(yY , =7r + foarccos h + L 



2 , -2 



(b + LY + ir 

At the same time, in Ref. |L2| it has been shown that the exact three-loop solution 
(with Pade transformed beta-function) can also be expressed in terms of the Lambert 
function 



n (<r\ - 1 1 r _ --lnx/6+c-l . _ AsA (k \ 

a3{x) -~b' l-c + W(z)> Z ~~ e ' C ~^f 



In what follows, referring to the three-loop case, we shall imply the MS scheme with 
Pf 5 = 2857/2 - 5033^/18 + 325n^/54 . Here, at n f = 3, we have (3 3 = 3863/6 = 
643.833 ; c= 1.415. 



Besides the pole at x = 1, it has a logarithmic branch point at i c = e . 
6 The 3-4-4 per cent difference can be "compensated" by — 8% correction of the A value. 
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2.3 Stability of analytic coupling in the IR region 

The analytic coupling A(x) obeys several important properties in the IR region. It is 
quite stable, at Q < A, with respect to higher-loop contribution and to renormalization 
scheme dependence. Its IR limiting value 

.4(0) = 1 (9) 

(corresponding to a Sian (0) = 47r//? 1 (3) ~ 1.4) is universal in this sense. 

This remarkable property of universality was first noticed in our starting up paper M . 



Later on, a detailed analysis [g, 14, 11, 12 has revealed that the -4(0) value turns out to be 
insensitive not only to higher loop terms in the beta-function but to the precise structure 
of the ghost singularity (removed by analytization) as well. This structure depends on 
approximation. 

For instance, instead of the square- root singularity of the two-loop exact solution ([5]), 
an iterative approximate solution Eq. @ obeys a pole at x = 1 and a log's branch point 
at = e~ b . At the three loop case, the Pade-approximated solution obeys a pole and 
a branch point (discussed in detail in the paper flr2|), to be compared with the (lnx) -1 / 3 



singularity of the solution with a non-transformed beta-function. 

Nevertheless, in all these cases, the final analytic results for A(x) obey the property 
([]) and their IR behavior in the interval (0, 1) is very close |TD[ to each other. The 
analytization procedure "smoothing over all sharp angles" makes all them equal. 

One can also say that the spectrality bounds from above the invariant coupling by 
this maximal value (0), "keeping it reasonably small". For instance, usual beta-function 
in the MS scheme for nj = 3 numerically, 

/3(a) = - a 2 (l + 0.790a + 0.883a 2 ) a 2 (l + 0.566a + 0.453a 2 ) . (10) 

looks quite good: its higher terms are reasonably decreasing in the usual physical region 
with a < 0.4 (a < 0.3). It is worthwhile here to introduce the "beta-function for the 
analytic coupling" 

which is a trancendental nonanalytic function (3 an {A) of its only argument (with the IR 
fixed point at A = 1). In the one-loop case this function can be analyzed rather simply 
- see, below, Eqs.(|14]) and (|15"D. It turns out that its maximum value /3 an (l/2) = 1/12 is 
of the same order of magnitude (this statement remains valid in the higher loop case) as 
the expression flT0|) taken at a = 0.3 . 

Let us remind to the reader, that the connection between analyticity of a function in 
the cut complex plane and its boundness have been discussed several decades ago (see, 
e.g., Refs.fi"5|) in the context of a low-energy behavior of hadron scattering amplitudes. 
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2.4 Analytization of observables 



In papers || ||| a specific recipe for analytization of an observable M(s) has been in- 
troduced. First, one should relate M(s) possibly given in the time-like region, to some 
auxiliary function f(Q 2 ) of a space-like argument^ which obeys the analyticity property 
in the Q 2 plane compatible with a spectral representation of the Kallen-Lehmann type. 
This function /, after usual RG machinery, acquires a form of perturbative power series 

F(a(x)) = Y,fna n (x). (12) 

n 

By prescription first used in paper ||, the causal analytization means 

F(a(x)) F(x) = J2fnMx) (13) 

n 

with A n (x) = [a"(x)] an being "the n-th power of a(x) analytized as a whole". Here, 
each A n (x) satisfies the Kallen-Lehmann representation with the spectral density p n (o") 
defined as Sa n (— a). 

Note that expansion (|T3|) is not a power one as A n (x) ^ [A(x)] n at n ^ 1. The recipe 
(P^) =>- ([13]) changes the nature of expansion ! 

The surprise feature of this particular recipe called the "Analytic Perturbation Theory" 
(APT), is the remarkable stability of its results with respect to a higher loop contribution 
H ^] and, in turn, with respect to a scheme dependence |3T| . We are going to show 



that the origin of these physically important properties lies in the change of the type of 
perturbation expansion. 

The non-power set {A n {x)} is an asymptotic one at UV as 

A n+ i(x) ~ [lnx]~ (n+1) = o(An{x)) , x -> oo . 

This means that the series (|13D, generally, should be treated as an asymptotic expansion 
of the function J^(x) over an asymptotic set {A n (x)}. Its convergence features are 
determined, on the one hand, by the coefficients f n calculated on the basis of n-loop 
Feynman diagrams and, on the other, by the property of the set {A n (x)}. Turn to the 
discussion of this set(s). 



3 Set of expansion functions 
3.1 One- loop expansion functions 

The simplest set { [a n (x)] an } = {*4. n (x)} consists of analytized powers of the one-loop 
pQCD expansion parameter a(x). Besides (0) it includes 

M^) = ^- P 2(x)- P2 ( X) =-^— =P2 (iy 1 (M) 

7 As, e.g., the e + e~ annihilation cross-section ratio R(s) is related to the Adler function D(Q 2 ). 
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A I \ 1 X I X .,. 1 n / \ P2(x) 
A 3 (X) = ^-g- + T rg - -z-r ^ \ A 4 (x) = — j p 2 (x) — 

In x (1 — x) A 2(1 — x) z in x 6 

These "analytized powers" obey a specific symmetry 



= (-l)M n (l/a:) ; (n > 1) 
and are related by recursion relation with the help of the operator D = —x(d/dx) 

A n+1 (x) = -VA n (x) = ±-V n+1 \nA (x) with A (x) = ^± (15) 
n n! xmx 

being the generating function. 
In particular, 



(16) 



= ee (a an (x)) = /?« (a an (l/x)) . 

Here, P2n( a ) * s a non-analytic function of its argument. 
An attempt to express A n (x) via A(x) = a an (x), the analytic coupling, gives 

A 2 (x) = A\x) - ^—A{x) + — !— ; (17) 
1 — x I — x 



A 3 (x) = A\x) - -^A\x) + 3 -4(x) - ^ + 2 
I — x 1-ir 2(1 — a; 



— a sort of a "mixed" representation, combining polynomial and nonperturbative (via 
x = Q 2 /A 2 argument) dependencies. It can be argued (see below Section |3.2| ) that this 
representation is not interesting from a pragmatic point of view. 

In the two-loop case, to relate solution (|6|) with higher analytized powers, one can 
formally use the operator 



^ 1 d 1 + W(x) d 

V-2 = TT^TT = rTTT—, X 



so that A { *Ux) = ^A^(x) . (18) 



bdW W(x) dx n+iV ' n\ 

Note also, that D in Eq.(|I5|) can be treated as an operator of differentiation over an 
"effective time variable" t = In x. An analogous interpretation of T> 2 gives 



„ 1 d d 
T> 2 = = — 

l + fea 2 (e*) dt dr 



r = t + b /*a(e*'W, (19) 
J o 



with r(t) = t 2 , an "effective two-loop time". For large t values one has r ~ t + feint. 
However, expression flT8| ) needs to be further specified at < x < 1. 
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3.2 Subtraction Structures and behavior at "the low Q region" 



The rational structures p n (x) = A n (x) — a n (x) that subtract ghost singularities 

Pl{x) = ih ; Mx) = (T^F ; P3{x) = w^w ; 

P*( X ) = - X{X ~ 6 ( 1 + 1%1 X ~ ) > (^ = 2±V3); ... (20) 

are connected by a recursion relation analogous to (|T5[) and, except p±, obey (anti) symme- 
try under x — > 1/x . 

As can be seen from this recursion relation, all p n >2{x) at the origin have the first-order 
zero that provides a property 

An(0) = ; n > 2 (21) 

valid in the higher loop case as well. 

For a quantitive orientation it is useful to study the A n {x) behavior at a; = 1+e : e < 1. 
At the one-loop case this can be done explicitly with the help of (|2~0"D 



A(x) ~ - - — ; A 2 (x) ~ — + — e 2 ; A 3 (x) ~ H — — ; ^ — — + — — 

K J 2 12' y 1 12 60 ' K J 240' V ; 720 3042 



These numerical results are rather instructive. Together with Eq. (|2lD they show that 
in the "low Q region", at Q < A, we have 

\A n (x)\ <£A n (x), (22) 

an important estimate, which, at the very end, is responsible for a low level of sensitivity 
of APT results with respect to the higher loop and scheme effects - see, e.g., observations 
made in Refs. || - fpTOfl . The last estimate is valid also in the two- and three-loop cases. 



It is not practical to use expressions like (||) , (§) and fllS]) for explicit analysis of asymp- 
totic sets {^4 n (a;)} at two- and three-loop cases. For a short quantitative discussion we 
rather use results of numerical calculation via an adequate spectral integral. They show 
that in the interval < Q < A in addition to relation (|22"D we have 



Ak+i(x) = [A{x)\ 



l+2fe 



In other words, numerically, in the "low Q region" the real expansion parameter is closer 
to A 2 (x) rather than to A(x) Q. Naturally, as Q/A grows and power terms diminish, all 
A n (x) tend to their natural limits |/4.(a;)] n ~ a n (x). 

Moreover, as it follows from the representation ([20) and property of singularity In - " - x at 
x = 0, the expansion functions A n+ 2(x) obey precisely n zeroes on the interval (0, X(A)) 

8 In particular, this means that a "mixed" representation in powers of A(x) with nonperturbative 
cocmcients is not reasonable due to big cancellation inside the r.h.s. of relations analogous to Eq.(p7|). 
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with X(A) being the upper boundary of the region where a power nonperturbative cor- 
rection ~ x -1 = A 2 /Q 2 is essential. Hence, the set under discussion consists of quasi- 
oscillating functions. This feature makes the problem of estimating the resudial term 
(that is an error) in the asymptotic expansion Eq. flT3]) more complicated. This is typical 
for the asymptotic expansion a la Erdelyi. 

Note also that the sets {^4 n (a;)} both in the one- and two- loop cases obey a peculiar 
structure. Their neightbouring terms are connected by differential relations fll5D and flT8|). 



4 Discussion 

We have analysed a particular version of "Causal, Q 2 -analytic perturbation theory" , the 
APT version, which, by convention first introduced in paper 0, uses a set {A n (x)} for 
analytization of observables. It can be considered as a "nonpower analytization" , to 
distinguish it from another possibility, the "power" one with the help of an asymptotic 
set {A n (x)} by the recipe 

F (o(x)) =► F (A(x)) = £ fn [A{x)f , (23) 

n 

instead of ([X3|) . 

Just the nonpower analytization yields intriguing results with respect to loop and 
scheme stability. At the same time, the power analytization, Eq. (|2"3"D , results in a moderate 
change of usual pQCD practice mainly in the IR region. 

This, technically simpler, second version has an advantage from a theoretical point of 
view related to the issue of Consistency of analytization with the RG structure - see below 



Section |4.2| . To clarify, let us make a comment on the structure of the RG algorithm and 



on "noncommutativity" of analytization with some of its elements. 
4.1 On ambiguity of analytization procedure 

The procedure of the renorm-group method, in addition to deriving functional and diffe- 
rential group equations, consists in a few steps 
[1] Calculating beta-function(s) and anomalous dimensions; 
[2] Solving RG differential equations (RGDEs) for invariant coupling(s) a(x); 
[3] Solving RGDEs for other functions f(Q 2 , a), e.g., propagator amplitudes, effective 
masses and "physical amplitudes''^ with the use of explicit expressions for beta-function 
of Step [1] or invariant coupling(s) a(x) obtained in Step [2]. The resulting F (a(x)) can be 
expressed as a power series (starting, possibly with logarithmic term). 

The invariant coupling analytization adds an additional step that follows the Step [2]: 
[2a] a(Q 2 )^A(Q 2 ). 

9 Like , Adler functions, structure function moments, etc. 
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However, analytization of propagators and observables can now be performed either by 
modification of Step [3] - 

[3m] Using explicit expression A(x) in the process of RGDEs for f(Q 2 , a) solving 

or as an additional Step: 
[4apt] Analytizing the result of Step [3], i.e., by applying the analytization procedure to the 
power series for F (a(x)) =>• F(x). 

The sequence 

[1] + [2] + [2a] + [3] + [4 APT ] = [APT] 

was used in Refs.f|-[|l(J. Just this procedure yields nonpower asymptotic expansion flllf ). 

On the other hand, in parallel with step [4apt] there exists a simpler possibility: 
[4 an ] Substituting expression A(Q 2 ), like in (p3|), in the result of Step [3]. 

The sequence 

[1] + [2] + [2a] + [3] + [4J = [ICA] 

can also be used for analytization of observables. This procedure, involving just the 
'invariant coupling analytization' (ICA), yields power asymptotic expansion fl23|) differing 
from the usual one, Eq.([L2D, by substitution a(Q 2 /A 2 ) A(Q 2 /A 2 ) only. 

We see that, generally, a causal analytization is not a unambiguous procedure. Quite 
remarkably, the above-mentioned ambiguity is of a functional nature (not in possibility 
to introduce an adjustable parameter). 

4.2 Analyticity vs RG structure ? 

Meanwhile, the sequence [1] + [2] + [2a] + [3m] contains an inner contradiction. E.g., 
Step [3m], used for the gluon propagator amplitude^ yields an expression that, at the 
very end, is not compatible with the result of the previous Step [2a]. At one loop level, 
it gives [|l7j 

<? RG {x) = [a(x)P =► [M*)P 

with Ao(x) defined in fll^P. However, as it follows from basic RG relations, the product of 
a vertex and appropriate powers of propagators forms an invariant coupling. In the case 
under consideration, one obtains Aq{x) rather than a an (x) = A(x) used as an input. 

Quite analogously, there is a subtlety with the Step [4apt] implementation. The point 
is that for some objects (e.g., for propagator amplitudes) the result of Step [3] at the one 
loop level starts with fractional power (or logarithm) of a(x) that gives rise to a branch 
point. The analytization of expression like [lnx] _I/ is equivalent to subtraction of a cut 
contribution, i.e., yields a two-term structure of a specific form[]]. It is easy to see that the 
appropriate product of such structures cannot give expression (H) for invariant coupling. 
Taken literally, this observation means that the APT procedure also faces a contradiction 
with the RG structure. 

10 Compare with the Step [4 an ] used in Ref.Q. 
11 For explicit expressions we refer, e.g., to Refs. || 
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Conclusion 

1. Our analysis in Section 3 reveals that the APT expansion, Eq.(^), for an observable 
function, generally, represents an asymptotic expansion over a nonpower asymptotic set 
{A n {x)}. The latter obeys quite different properties in various ranges of the x variable. 
In UV it is close to the power set {a n (x)}, commonly used in the current practice of QFT 
pertirbation calculation. Hence, the APT converegence property in UV is completely 
determined by expansion coefficients f n . On the other hand, in IR the asymptotic set 
{^4 n (x)} is of a more complicated structure. In the "low Q region" the behavior of the 
functions A n (x) ;n > 3 is oscillating. Due to this, the contribution of higher terms in 
the APT expansion is suppressed. The APT expansion, Eq . (|T3"D , in IR has a feature 
of asymptotic expansion a la Erdelyi. This tentative conclusion raises hopes that the 
pertubative approach to QCD may be fruitful in the region Q ~ 1 GeV where the QCD 
running coupling is not a small quantity. 

2. In Section |4.1| , we have shown that the general program of Invariant Analytisation, 
being quite a definite procedure for effective coupling, is not "rigid" enough when applied 
to other objects. In particular, it contains a degree of freedom in analytizing observables. 

This ambiguity together with a "proximity to contradiction" , discussed in Section |4.2j . 
poses a question of looking for an additional ansatz in the whole Causal analytic approach. 
The APT possibility is too interesting to be "abandoned without a struggle" . 

3. In our opinion, one more funny lesson of the considered nonpower construction is 
a semiquantitive observation that the APT approach is equivalent to the usual pQCD 
practice with one strange amendment: "To restrict calculation to only the leading QCD 
contribution" ; by the way, forgetting about all headaches of higher-loop diagram calcu- 
lation, scheme dependency and expansion convergence. This intriguing feature could be 
formulated as a suspicion that the pQCD is an effective theory^ and its higher order 
contributions have no clear physical content. 
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